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Abstract 

Applying a modified version of the Gledzer-Ohkitani-Yamada (GOY) shell model, the signa- 
tures of so-called two-dimensionalization effect of three-dimensional incompressible, homogeneous, 
isotropic fully developed unforced turbulence have been studied and reproduced. Within the frame- 
work of shell models we have obtained the following results: (i) progressive steepening of the energy 
spectrum with increased strength of the rotation, and, (ii) depletion in the energy flux of the forward 
forward cascade, sometimes leading to an inverse cascade. The presence of extended self-similarity 
and self-similar PDFs for longitudinal velocity differences are also presented for the rotating 3D 
turbulence case. 
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I. INTRODUCTION 



Rotation of a fluid has been used as a control parameter that can progressively make 
a 3D turbulent flow look like a quasi-2D or a 2D turbulent flow. This phenomenon - 
known as two-dimensionalization of three-dimensional (3D) turbulence — may be termed 
two-dimensionalization effect. The phrase 'look like' in this connection basically means that 
certain properties of 3D turbulence, such as wavenumber dependence of energy spectrum, 
direction of energy cascade etc., change to such behaviors that are more recognized with a 2D 
turbulence flow; in other words the flow would seem to have become 'two-dimensionalised'. 
In view of the fact that the dynamics of oceans, atmospheres, liquid planetary cores, fluid 
envelopes of stars and other bodies of astrophysical and geophysical interest do require an 
understanding of inherent properties of turbulence in the rotating frame of reference, the 
problem of two-dimensionalization is naturally of profound scientific importance; turbulence 
in rotating bodies is even of some industrial and engineering interest. 

Two-dimensionalization appears to be a subtle non-linear effect, which is distinctly different 

n n 

from Taylor-Proudman effect, as shown in the works of Cambon[l|], Waleffe[2[ and others. 
Numerical simulations js]] indicate the initiation of an inverse cascade of energy with rapid 
rotation, a fact well supported by 
Baroud et al (4, 6] and Morize et al. 



;he experiments |4j, [5] . Although recent experiments by 
7| have shed some light on the two-dimensionalization 



effect, the scaling of two-point statistics and energy spectrum in rotating turbulence remains 
a controversial topic. An energy spectrum E(k) ~ k~ 2 has been proposed js, 9] for rapidly 
rotating 3D turbulent fluid and this does seem to be validated by some experiments 4], Q| 



and numerical simulations 



10 



]Hl3[. But some experiments [5|] do not agree with this proposed 
spectrum. They predict steeper than k~ 2 spectrum and this again appears to be suppor ted 



from by results 



14 



15] and analytical results found using wave turbulence theorv|16l. Il7| 



n 



It is well known that structure functions for 3D [18], quasi-2D[19| and 2D 
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211 turbu- 



lences contain quite a lot of information about the respective flows; for instance the exact 



results for t 



lird-order structure functions serve as benchmarks for any theory of turbulence. 



Recently 22j, |23| , there has been an attempt to calculate such non-trivial results for rotating 
turbulent flows. Based on those results, it has been argued |24j that the presence of helicity 
cascade in the rotating flow would cause depletion in the forward cascade of energy that 
sometimes may lead to inverse cascade and that the exponent (— m) in the energy spectrum 



E{k) ~ k m should lie between —2 to —7/3 for rapid rotation. 



In this paper, we shall use GOY shell model[25|, |26j, modified appropriately, to investigate 



these signatures of two-dimensionalization effect, the behaviour of the structure function 



and the status of extended self-similarity (ESS)[27J in the rotating flows 



One may indeed ask why one needs another shell model after Hattori et. al. [11] already have 
proposed modified version of Sabra shell model 28[ a few years ago. To answer this question, 
let us collect the main results of that model: i) the exponent (— m) of the energy spectrum 
in the inertial range changes from —5/3 to —2, ii) no inverse cascade is detected with the 
increase in rotation rate, and iii) the PDF's of the longitudinal velocity difference doesn't 
match with the experiments. Studies of the last few years' research-literature on turbulence 
would reveal that investigations of two-dimensionalization effect are growing rapidly. Neat 
experiments have confirmed that the exponent (— m) overshoots the value —2 quite com- 



fortably as for instance shown in the experiments by Morize et. al. [5j. Moreover, the fact 
that, some experiments and numerical simulations do show inverse cascade with increase in 
the rotation rate, motivates us to construct shell models that can mimick this effect. As 
mentioned above, Hattori et. al.'s model finds PDF which mismatches with experiments 
and also, the model requires a fluctuating part in the rotation rate to obtain various results, 
while in experiments and simulations such effects are not known. This again should motivate 
the need for another model. Moreover, the numerical simulations performed in the present 
paper are mainly for decaying turbulence while the model of Hattori et. a/.'s dealt with 
forced turbulence. Hence, in this paper we try to consider another possible shell model that 
can mimic the signatures of the two-dimensionalization effect in details. 

II. SHELL MODELS: A BRIEF OVERVIEW 

Shell models for fluids are, in practice, simplified representative versions of the Navier- 
Stokes (NS) equations; but they do retain enough of the flavour of the parent equations 
making themselves handy testing grounds for many statistical properties of fluid turbu- 
lence. In fact, shell models have been used to study statistical properties of turbulence 



in the past 



29 



301 ] with a fair degree of success. The fact that 3D turbulence still lacks 
solid understanding, can be related to the lack of complete theoretical characterization and 
explanation of the energy- cascade mechanism — the process that spreads and sustains tur- 



bulence over wide range of scales. The popularity of shell models is due to their usefulness 
in modeling this very energy- cascade mechanism. The other advantage of using shell models 
is that, being deterministic dynamical models, they can be studied by dint of faster and 
accurate numerical simulations. This stems from the fact that the number of degrees of 
freedom (DOF) needed to reach high Reynolds numbers(-Re) is just moderately high as the 
number of DOF grows logarithmically in Re as opposed to NS case, where number of DOF 
goes as Re 9 / 4 . 

In brief, shell models are a set of coupled nonlinear differential equations each labeled by 
index 72 = 0,1,2,..., called the shell index: 

+ vk 2 ^j u n = k n (NL) n [u, u] + f n (1) 

where the complex dynamical variables u n represent the temporal evolution of velocity fluc- 
tuations over a wavelength k n . The wavenumbers k n are given by k n = ko\ n with A, called 
the intershell ratio, usually set to 2 and ko being a reference wavenumber. The forcing 
term is taken as time-independent and is usually restricted to a single shell: f n = fS nn *. 
Velocity evolution is thus followed over a set of logarithmically equispaced shells. The non- 
linear term, (NL) n [u,u\, is so chosen that in some sense total energy, helicity and phase 
space volume are conserved as is typical of the nonlinear term in the NS equation. Another 
usual practice is demanding locality of interactions (only between nearest and next nearest 
shells) in shell(Fourier) space, although this is not absolutely necessary. The advantage in 
doing so is one gets rid of the so called sweeping effect (the direct coupling between inertial 
and integral scales) making shell models ideal to study nontrivial temporal properties of the 
energy cascade mechanism because the time fluctuations are not obscured by the large-scale 
sweeping effect. This makes shell models rather approximate, quasi-Lagrangian representa- 
tions of NS equations. 

The choice of (NL) n is not at all unique and thus quite a few shell models have been 
proposed in the literature. Most popular of the lot being the so-called Gledzer-Ohkitani- 



Yamada (GOY) model 25|, |26(]; another example being the sabra model 281 ] . We however will 
concentrate solely on the GOY model and modify it appropriately to incorporate the effect 
of rotation. 



III. THE MODEL 



The equation of motion describing 3D turbulence are the Navier-Stokes equations. For 
an incompressible fluid in a rotating frame it reads as: 



— + u • Vu + 2S1 x u 

at 

V- u 



--Vp + z/V 2 u + f 

P 







(2) 
(3) 



where O is the angular velocity of system rotation, v the kinematic viscosity, p the density 
and f the external force. The term 2f2 x u in the equation is due to the Coriolis force 
and thus is absent in the non-rotating case. As an aside, it may be mentioned that such 
type of linear terms can also be seen in complex cascade models for magnetohydrodynamic 



turbulence 



ail. 



In this paper, we have adopted the following strategy for the numerical investigations 
A specific form of GOY shell model for non-rotating decaying 3D turbulence is: 
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33). 



d 1 2 



u, 



ikr, 



U n +2U n +i — —U n+ iU n -i — — U n _i« n _2 



(4) 



As discussed before, this may be thought as a time evolution equation for complex scalar 
shell velocities u n (k n ) that depends on k n — the scalar wavevectors labeling a logarithmic 

We choose: ho — 1/6, v — and n = 1 to 22. The 



discretised Fourier space (k n = k 2 



initial condition imposed is: u n = kn' l e ldn for n — 1, 2 and u n = kn 2 e~ h ™e tBn for n = 3 to 22 
where 6 n G [0, 2ti} is a random phase angle. The boundary conditions are: u n = for 
n < 1 and n > 22. It is easy to notice that GOY model has certain similarities with the 
NS equation. The nonlinear terms have dimension [velocity] 2 / [length] as is the case with 
the nonlinear term(u • Vu) in the NS equation. In the inviscid limit {y — > 0), equation (jl]) 
owns two conserved quantities viz., J2 n \ u n\ 2 and J2 n (~ ^) n ^n\u n \ 2 which are seen on the 
same footing as energy and helicity respectively. Apart from this and most importantly, this 
model displays multiscaling, i. e. it has been shown 3^] that at inertial scales, the structure 
functions (may be naively defined for shell models as (\u n \ p )) display power law dependence 
on k n with non-trivial exponents: {\u n \ p ) oc k n ^ p , where the exponents ( p have a non-trivial 
nonlinear dependence on the order p. 

If the fluid is rotating then one may modify equation (jl]) by adding a term R n = 
— i [u + (— l) n h]u n in the R.H.S. This specific form of the term R n was originally intro- 
duced by Reshetnyak and Steffen[12j. u and h are real numbers. It may be noted that 



FIG. 1: A representative plot of structure functions S p vs. k n , on log-log scale for uj = 0.01, h = 0.1 
for E p vs. k n . From the topmost curve to the bottommost curve p increases from 1 to 6. We plot 
for n = 3 to 20. 

this term, as is customary of Coriolis force term in NS equations, wouldn't add up to the 
energy. The (— l) n h part in R n has been introduced to have non-zero mean level of helicity 
that otherwise has a stochastic temporal behaviour and zero mean level. Therefore, the 
appropriate shell model for rotating 3D turbulent fluid is: 

1 1 



d r 2 




It + 





U n +2U n +l — ~;U n+ iU n -i — -U n -\U n -2 
4 o 



-i[oj+(-l) n h]u n (5) 



We fix h = 0.1 in our numerical experiments and test for the following values: u> = 
0.01,0.1,1.0 and 10.0. Later, we shall come back to the effect of changing h while keep- 
ing uj fixed. We found ourselves poor in computational resources when trying to simulate 
for higher u, say 100. the practical problem is that higher the u, longer one has to run the 
simulation to get an appreciable — say, 10-shell-wide — inertial range. We shall henceforth 
refer u as rotation strength. 

All the data points reported here are averaged over 2000 independent initial conditions and 
the error-bars reported herein are the corresponding standard deviations. Data have been 
recorded only after the energy cascade has stabilized and a nice the inertial range can be 
comfortably defined between shell n = 4 to 15. We have applied the slaved second order 



Adam-Bashforth scheme 



351 ] to numerically integrate equations @ and ([5]). 



The pth order equal time structure function (see Fig. 1) for the model has been defined as: 

1 



Im 



(6) 



to avoid period three oscillations 36|. The energy spectrum has been defined as: 

E(k n ) = Y, p {k n )/k n ~ k- m . (7) 
The mean rate of dissipation of energy is, of course, 

(8) 

and flux through nth shell is calculated using the relation: 

Il n = ( -Im k n u n+1 u n (u n+2 + \u n -i J ) (10) 



e = l^^vk 2 n \u n \ 2 ^j 



4 

For studying relative structure function scaling, the ESS scaling exponents are taken as 

c; = c P /c 3 . (ii) 

The exponents m, ( p and Q have all been estimated for inertial ranges only 

IV. THE RESULTS 

In this section, we systematically present the results obtained by the numerical simula- 
tions and based on these results we show, that indeed the two-dimensionalization effect is 
in fact mimicked by the present shell model. 

A. Signatures of two-dimensionalization 

One can clearly see, Fig. 2 and Fig. 3, that as the rotation strength increases, the energy 
spectrum becomes steeper and the slope monotonically increases from a value ~ —5/3 to a 
value of ~ —7/3; hence validating one of the two-dimensionalization effect's signatures. 
Investigating the direction of the flux in the inertial range regime, we find (see Fig.4) that 
with the increase in rotation strength, then first the forward cascade rate starts decreasing 
and furthermore instances appear when at certain shells the flux direction reverses. Again, 
the number of shells with such behavior increases as the rotation strength is enhanced; 
clearly suggesting a depletion in the rate of forward cascade. Thus, yet another signature of 



Table 1: ( p for p = 1 to 6 for various rotation strengths. 


P 


(p(lj = 0.00, h = 0.0) 


C p (u} = 0.01, h = 0.1) 


( p (u = 0.10, h = 0.1) 


C P (w = 1.00, h = 0.1) 


Cp(w = 10.0, h = 0.1) 


1 


0.37 ± 0.0027 


0.52 ± 0.0086 


0.63 ± 0.0098 




0.62 ± 0.0067 


0.66 ± 0.0086 


2 


0.70 ± 0.0062 


0.95± 0.0182 


1.1 ± 0.0232 




1.2 ± 0.0161 


1.2 ± 0.0138 


3 


1.0 ± 0.0127 


1.3 ± 0.0394 


1.6 ± 0.0455 




1.7 ± 0.0301 


1.8 ± 0.0197 


4 


1.3 ± 0.0251 


1.7 ± 0.0712 


2.0 ± 0.0733 




2.2 ± 0.0490 


2.4 ± 0.0283 


5 


1.5 ± 0.0454 


2.0 ± 0.1083 


2.3 ± 0.1017 




2.7 ± 0.0713 


2.9 ± 0.0402 


6 


1.8 ± 0.0718 


2.4 ± 0.1470 


2.7 ± 0.1291 




3.2 ± 0.0953 


3.4 ± 0.0550 


Table 2: £* = Cp/Ci f° r p = 1 to 6 for various rotation strengths. 


P 


Cp(w = 0.00, h = 0.0) 


Q(u = 0.01, h = 0.1) 


Cp(w = 0.10, h = 0.1) 


Cp*( 


oj = 1.00, h = 0.1) 


Q(oj = 10.0, h = 0.1) 


1 


0.37 ±0.0153 


0.40 ± 0.0480 


0.39 ± 0.0553 




0.36 ± 0.0368 


0.37 ±0.0283 


2 


0.70 ±0.0188 


0.73 ± 0.0576 


0.69 ± 0.0687 




0.70 ± 0.0463 


0.67 ±0.0335 


3 


1.0 ±0.0253 


1.0 ± 0.0789 


1.0 ±0.0910 




1.0 ± 0.0603 


1.0 ± 0.0393 


4 


1.3 ±0.0377 


1.3 ±0.1106 


1.2 ±0.1188 




1.3 ± 0.0791 


1.3 ± 0.0479 


5 


1.5 ±0.0580 


1.5 ±0.1477 


1.4 ±0.1472 




1.6 ±0.1014 


1.6 ± 0.0598 


6 


1.8 ± 0.0844 


1.8 ±0.1865 


1.7 ±0.1746 




1.9 ±0.1255 


1.9 ± 0.0746 



two dimensionalization has appeared the shell model studies. 

At this point, it must be appreciated how important the inclusion of term — £(— l) n h in 
equation (jSJ) is in responsible for the effect of depletion in the rate of forward cascade. By 
setting mean level of helicity above zero, it is this very term that — in accordance with the 
arguments 23| that it is the helicity that is causing this signature of two dimensionalization 



effect to show up — has empowered the model with the capacity to mimic the effect: at- 
tempts to see this very effect when setting h = fails. To illustrate this, we fixed u at 0.1 and 
increased h monotonically. As Fig. 5 shows, for h = one doesn't get instances of negative 
flux; however, as h increase, several of such instances can be observed. We feel that for our 
simulations h — 0.1 is a good choice to arrive at various features of two-dimensionalization 
effect using this modified GOY shell model. 



FIG. 2: Energy spectra E(k n ) vs. k n plotted in log-log plot. Asterisk is the marker for non- 
rotating case whereas square, triangle, circle and diamond respectively are the markers for to = 0.01, 
to = 0.1, to = 1.0 and to = 10.0 cases. We plot for n = 3 to 20. For clarity, in an accompanying 
figure, we have also plotted compensated energy spectra k^E(k n ) vs. k n only for non-rotating and 
to = 10.0 cases. One may note how the slope changes from 1/3 to —1/3 with rotation as has been 



predicted 
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FIG. 3: The slopes of the energy spectra (obtained from Figj2|) plotted against the "so-called" 
rotation strength. The accompanying dashed line is the value —7/3 of the slope that has been 
predicted for very rapid rotation. 

B. Extended Self-Similarity 

The study of Extended Self-Similarity (ESS) in the shell model has also been revealing. As 
can be seen in Fig. 6 and Tables 1 and 2, the increase in the rotation strength is accompanied 
by a departure from the usual She-Leveque scaling. But, the fact that at higher p, ( p 



4 5 6 7 8 9 10 11 12 13 14 15 
n 



FIG. 4: Average flux of energy through nth shell vs. shell number n. The function Tl n at the 
y-axis helps to plot all the curves distinctly in a single figure. Only the inertial range (n = 4 to 
15) has been plotted. Markers are same as that for Fig. 2. 




4 6 8 10 12 14 



n 



FIG. 5: Average flux of energy through nth shell vs. shell number n. Only the inertial range 
(n = 4 to 15) has been plotted. Asterisk, square, triangle and circle respectively are the markers 
for h = 0.0, h = 0.01, h = 0.1 and h = 1.0 cases, to has been kept fixed at 0.1. 

seemingly becomes parallel to p/2 vs. p, is worth paying attention: This is in accordance 
with the direct numerical simulation (DNS) resultspj^ and experimental results^]. However, 
most interestingly is probably the observation that within the statistical error, Q = Cp/Ca 
obtained for the rotating system via ESS coincides with that for the non-rotating ones. 
Probably, this extends the ESS for 3D fluids even further by implying that rotation keeps ESS 
scaling intact, even though usual ( p changes owing to rotation. Of course, only experiments 
and DNS can judge if this really is true for real fluid turbulence: The GOY shell, is after all 
just a model that remarkably well reproduces many characteristic features of turbulence by 
only using a fraction of computation power needed by DNS. In this context, one might be 
well aware that some modified versions of GOY model invented to describe the distinguishing 



1 2 p 3 4 5 6 



FIG. 6: Cp vs - P plotted for the data in Table 1. Markers are same as that for Fig. -2. The dashed, the 
chain and the dotted lines are respectively for ( p = p/3 (K41), ( p = p/2 and ( p = p/9+2[l-(2/3) p / 3 ] 
(She-Leveque exponent \s7\). The dotted curve has almost been reproduced by non-rotating GOY 
model, as expected. This anomalous scaling is remarkably reproduced in the model dynamical 
system with limited number of degrees of freedom because its chaotic evolution exhibits temporal 
intermittency[34]. The inset is plot for £* vs. p plotted using the data of Table 2. All the connecting 
lines and the fractional values of p are just aids for the eyes. 

n 

features of 2D turbulence have been shown to give spurious results |38|. Thus, one always 
has to be careful while dealing with simplified models of turbulence. 



C. Probability Distribution Functions 



We have also tried to see if we can get self-similar probability distribution function (PDF) 
for longitudinal velocity differences as has been reported in experiments^]. The GOY model 
is defined in /c-space but we study the aforementioned PDF in real space obtained by using 



a sort of inverse Fourier transform 



of the form: 



N 



v(r 



t) = ^c n [u n (t)e 

71=1 



+ c. c 



(12) 
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FIG. 7: (Color online) Normalized probability vs. normalized longitudinal velocity difference for 



been taken to be 0, 1, . . . , 7 and they respectively correspond to red circles, green asterisks, yellow 
squares, cyan stars, magenta diamonds, blue crosses, black triangles and orange pluses. 

Here, the wavevectors are k n = k n e n where e n is a unit vector in a random direction, for 
each shell n and c n are unit vectors in random directions. We ensure that the velocity field is 
incompressible, V -v = 0, by constraining c n -e n = 0, Vn. In our numerical computations we 
consider the vectors c n and e n quenched in time but averaged over many different realizations 
of these. 

Thus, when in the simulations energy cascade is stabilized, we use shell-velocities to find 
the real-space velocities following the aforementioned prescription. From the velocity field 
obtained in this manner, one can easily construct longitudinal velocity difference as: 



by dividing Sv\\ by their rms value. In Fig. [7j we present PDFs of normalized 5v\\ for 
decaying rotating turbulence with an cu-value equal to 10. One may note that the PDFs 
are non-Gaussian but, when re-scaled appropriately, are fairly self-similar — the plots for 
various separations (1) collapse on a single curve. 





(13) 





FIG. 8: (Color online) Normalized probability vs. normalized longitudinal velocity difference for 
the forced case. Markers follow the convention used in figure ([7]). 

However to make contact with experiments^], we must also study a forced version of the 
model. Hence, we repeated our simulations for the case of rotating GOY shell model (given 
by equation (JS}) with a forcing term in the R.H.S: 



d 1 1 



ik n 



U n -lU n _ 2 



-i[cj+(-l) n h]u n + fd n ^U) 



For this particular case, we chose / = 5 x 10 _3 (1 + i) and we increased the total number 
of shells from 22 to 24. The corresponding PDFs for u = 10, as given in Fig. (JSJ), are 
non-Gaussian but also nicely self-similar. Thus, the rotating GOY model can reproduce 
even the self-similar feature of the PDFs quite impressively. 

V. DISCUSSION AND CONCLUSION 



Shell models have been successfully used to study statistical properties of turbulence by 
many authors (see Ref. ( 29]) and Ref. ([3]) for details). Most of the studies have dealt with 
the case of homogeneous and isotropic turbulence. Hattori et. al. proposed a shell model 
for rotating turbulence. Here we have attempted to improve their results by investigating 
two-dimensionalization effect by using a modified version of GOY shell model. Some results 
of the model are, no doubt, consistent with experiments and DNS. 

Concerning our main aim — modeling the two-dimensionalization effect — one can always 
question the robustness of the obtained signatures because i) a scaling law for a single- 



component spectrum, though heavily used in literature, has poor meaning in the strongly 
anisotropic configuration which is relevant when passing from 3D-2D; different power laws 
can be found in terms of k z , k± and k in contrast to the 3D isotropic case, and ii) the inertial 
wave-turbulence theory is not consistent with an inverse cascade. Actually in weak-wave 
turbulence, getting rid provisionally of helicity and polarization spectra, a two-component 
energy spectrum e(k, cos 6) with cos6 = k z j a/ {k\ + k\) is found to be useful; if E denotes 
the traditional spherically averaged spectrum, the anisotropic structure is one of the best 
ways to quantify all intermediate states from isotropic 3D (with e = E(k)/(47rk 2 )) to 2D 
state (with e = E(k±)/(27rk±)8(k z )). Two-dimensional trends are therefore linked to a 
preferred concentration of spectral energy towards the transverse wave-plane k z = 0. This 



concentration, however, does not necessarily yield an inverse cascade 16|. A reasonable 
suggestion, in the light of this discussion, would be that in the shell model for rotating 
turbulence k should be interpreted as k±. 

In the closing, it may be concluded that this study has put the equation (jSJ) as a very good 
shell model for the rotating 3D turbulent flows; after all, it explains the observed signatures 
of the two-dimensionalization effect closely. Probably, this model and the model due to 
Hattori et. al. can together model the rotating turbulence in a simple but effective manner. 
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